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ABSTRACT 

The correlation function of the distribution of matter in the universe shows, at large scales, baryon 
acoustic oscillations, which were imprinted prior to recombination. This feature was first detected in 
the correlation function of the luminous red galaxies (LRG) of the Sloan Digital Sky Survey (SDSS). 
The final release (DR7) of the SDSS has been recently made available, and the useful volume is 
about two times bigger than in the old sample. We present here, for the first time, the redshift space 
correlation function of this sample at large scales together with that for one shallower, but denser 
volume-limited subsample drawn from the 2dF redshift survey. We test the reliability of the detection 
of the acoustic peak at about 100 /i~^Mpc and the behaviour of the correlation function at larger 
scales by means of careful estimation of errors. We confirm the presence of the peak in the latest data 
although broader than in previous detections. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: distances and 
redshifts — galaxies: statistics — large-scale structure of universe 



I. INTRODUCTION 

About 380,000 years after the Big Bang, when the 
temperature falls low enough for recombination to oc- 
cur, the matter in the Universe becomes neutral. At this 
epoch, the sound speed drops off abruptly and acoustic 
oscillations in the fluid become frozen. Their signature 
can be detected in both the Cosmic Microwave Back- 
ground (CMB) radiation and the large-scale distribution 
of galaxies. The characteristic length scale of the acous- 
tic oscillations can be used as a "standard ruler" to probe 
the expansion history of the universe. 

The imprints in the matter distribution of this feature, 
called baryon acoustic oscillations (BAO), should be de- 
tectable in both the correlation function and the matter 
power spectrum. Moreover, this feature should mani- 
fest itself as a single peak in the correlation function at 
about 100 /i^^MpcQ. The first unambiguous detection of 
this feature in the correlation function (redshift space) 
of the Sloan Digital Sky Survey Luminous Red Galaxie s 
(SDSS-LRG) was reported by lEisenstein et all (|2005l ). 
Detection of these oscillations in the power s pectrum of 
, the ga laxy distribution was first reported bv iCole et al.l 
(|2005D for the 2-degree Field Galaxy Redshift Survey 
(2dFGRS). Cosmological parameters can be determined 
from the position of the baryon acoustic peak (see, e.g., 
I^rcival et al. 2007). 

In this letter we study the correlation functions of the 
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largest available redshift surveys (2dFGRS and SDSS). 
We focus on the two-point correlation function, and 
study the reliability of the detection of the BAO fea- 
ture at 100 /i~^Mpc in the samples drawn from those 
surveys. In our calculations, we adopted a cosmological 
model with VLm = 0.27, and rJA = 0.73. 

2. DATA 

In the first d e tectio n of the baryon acoustic peak, 
lEisenstein et al.l (|2005l ) used a sample that was con- 
structed selecting about 15 luminous red (early- type) 
galaxies per square degree, u sing different luminosity 
cuts (see lEisenstein et "alll200lD . For correlation studies, 
an almost volume-limited (constant-density) subsample 
was chosen, its characteristics are listed in Table [1] as 
DR3-LRG. 

The last data release (DR7, see lAbazaiian et"ani2008| ) 
of the SDSS-LRG contains spectra for 206,797 LRGs 
within a solid angle of 9,380 square degrees. We se- 
lect a subsample of these data (the main compact body 
of the sky footprint, 7,204 square degrees) in order to 
minimize the influence of border corrections on our re- 
su lts. Choosing the same redshift and magnitude limits 
as lEisenstein et all (|2005l ) and slightly shifted magnitude 
limits (see Table 1) we have a sample that is about twice 
as large as the one used by them; we label it as DR7- 
LRG in Table [TJ This sample is approximately volume- 
limited, but not in a formal way. We hav e determined 
a rea l volume-limited (VL) sample (as in iZehavi et al.l 
120051 ). with nearly constant comoving number density, 
applying the appropriate luminosity cut and restricting 
the distance range to ensure completeness; it is labelled 
as DR7-LRG-VL. 

For comparison, we use a nearly volume-limited sam- 
ple from the full 2 dFGRS prepared by the 2dF team 
(jCroton et al.l 12004( 1 . It contains luminous galaxies in 
two spatial slices. The number density in the 2dFVL 
sample is more than an order of magnitude larger than 
the number density of the SDSS-LRG samples used here 
(see Table [TJ. 
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TABLE 1 

Characteristics of the samples used and quoted 



Sample Number Magnitude limits^ Redshift limits Solid Volume number density 

of galaxies angle (sr) (h^^Gpc^) (h^ Mpc^^) 



DR3-LRG 46,748 -23.2 < Mg(z = 0.3) < -21.2 0.16 < z < 0.47 1.16 0.75 6.3 x 10"^ 

DR7-LRG 102,568 -22.344 < Mg{z = 0.3) < -20.344 0.16 < z < 0.47 2.19 1.41 7.3 x 10^^ 

DR7-LRG-VL 44,164 -23.544 < Mg (z = 0.3) < -21.644 0.14 < z < 0.42 2.19 1.03 4.3 x 10"^ 

2dFVL 33,878 -21 < Aft, ^ < -20 0.03 < z < 0.19 0.45 0.023 1.5 x 10"^ 



Not e. — "Absolute magnitudes M are normalized to Hg = 100 km s ^ Mpc ^. ^ Using Qm = 0.3, Q\ = 0.7 as in lEisenstein et al.l 
II2005I '). these values are 0.72 and 6.5 X 10~^, respectively. 



Although the volume covered by the 2dFVL sample is 
smaller than that of the DR7-LRG, it is still useful to 
measure correlation on scales ~ 100 /i~^Mpc. As shown 
below, cosmic variance is not dominant, even if its effect 
is an order of magnitude larger for 2dFVL than for DR7- 
LRG. Moreover, the increase in density compensates the 
decrease in volume, so that the number of pairs of galax- 
ies in each distance bin at these scales is similar in both 
samples. Hence, discreteness errors should be similar in 
both cases. 

These three samples are important for the detection 
of the acoustic peak. Although the SDSS-LRG samples 
are larger and cover a redshift range less affected by non- 
linear effects, the lower redshifts mapped by the 2dFVL 
sample are also important, providing the yardstick that 
can be compared with the characteristic BAO scales at 
larger redshifts (see Fig. [T]). 

3. ESTIMATING CORRELATION FUNCTIONS 

We estimated the spherically-averaged redshift-space 
correlation functions by using the Landy-Sza lay border- 
corrected estimator (jLandv fc Szalavl ll993D that has 
good large-scale properties. We generated a random dis- 
tribution of points following the selection function of each 
catalogue considered, and estimated the correlation func- 
tion ^(s). 



Cls(s) = 1 + 



DD(s) DRis) 



RRis) RR{s) 



(1) 



where DD{s), RR{s) and DR{s) are the probability 
densities of galaxy-galaxy, random-random and galaxy- 
random pairs, respectively, for a pair distance s. There 
are several recipes for the choice of the size A^^d of the 
random point set; as we are interested in large-distance 
correlations, where the numbers of pairs per bin (kernel 
width) are large, we used N^d ~ {N is the number of 
galaxies in the sample) . Increasing N^^ up to 20A^ led to 
pointwise differences less than a percent. 

We estimate the probability densities by the ker- 
nel method, sum ming the box spline B^^-) kernels 
(|Saar et al.1 12007[) centered at each pair distance, and 
sampling the distributions at smaller intervals than the 
kernel width. 

For the 2dFVL sample, we generated the random cat- 
alogues using the subroutines to calculate the complete- 
ness and the magnitude limit fr om the angular ma sks, 
both provided by the 2dF team (jCoUess et al.ll2001[ ). 

For the DR7-LRG sa mples we gen erated the angular 
mask from the data (as inlHiitsiP2006l), following the scan 
stripes of the survey and defining by hand the rectangles 



in the survey coordinate system 77, A that cover the data 
footprint. We assumed the angular selection to be con- 
stant inside this mask. 

The DR7-LRG sample is only approximately volume- 
limited, and we generated the random samples based on 
the smoothed comoving density of the data. The ac- 
cepted statistical paradigm for the galaxy distribution is 
to model it as a Cox process - a Poisson point process 
where the local intensity is determined by a realization 
of a random field. The two-stage nature of this process 
leads to two sources of errors of sample statistics (cor- 
relation functions, in our case). The first source is the 
possible deviation of the realization correlation function 
from the true one (cosmic variance) and the other source 
is the discreteness of the point process - how well its 
correlation function estimate approximates that of the 
particular realization. 

Precise est imates of the cosmolo gical statistics have be- 
come actual; iNorberg et al.l (2008') classifies the error es- 
timates as internal and external, this coincides with our 
classification. As cosmic variance and discreteness vari- 
ance are independent, they add for the total variance. 

If the random field is Gaussian, as usually assumed, 
and its expected correlation function and power spec- 
trum are more or less known, the covariance of the cor- 
relation function estimate can be found as a convolution 
of the correlation function itself, or a s an integr al over 
the power spectrum squared (see, e.g.. lCohnll2006[ ). This 
is common for all samples; a rough upper limit is 



Var(^(r)) < 



P^{k)dk 



1 1 

2^VV2. 

5 • lO-^{V/h~^G-pc^)-\r/100h-^Mpc)-^, 



where the nume ri cal va lue was obtained by using the 
lEisenstein fc Hul (|1998f ) approximation for the power 
spectrum, with a Gaussian cutoff at fc = 100/iMpc~^ 
(the real-space scale of about 60/i~^Mpc). The sample 
volume V and pair distance r in this formula are typical 
for the DR7-LRG and the baryonic peak; the rms error 
is about 2 • 10^^ in this case, and about 1.5 • 10"'^ for 
the 2dFVL sample. As we shall see below, this is much 
smaller than the discreteness error in both cases, and we 
shall neglect it for the rest of the paper. 

The discreteness error can be estimated in several 
ways. The most attractive of these is bootstrap, that uses 
only the knowledge inherent in the observed data; the ob- 
served sample is repeatedly resampled and the statistics 
found averaging over these bootstrap samples. The as- 
sumption that the observed values are i.i.d. (independent 
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Fig. 1. — The large slice is drawn from the SDSS-LRG (DR7) survey. The slice is 6° wide in declination and the galaxy distribution is 
shown within the redshift range 0.16 < z < OAT. There are 10,136 red luminous galaxies within this slice depicted as red dots. The smaller 
slice with blue dots shows the galaxy distribution of 9,744 objects from the Southern Galactic hemisphere of the 2dFVL sample, reaching 
a depth oi z = 0.19. To illustrate the scale of the acoustic peak a segment of length 105 h^^Wlpc is shown to scale. 



and identically distributed) demands that resampling is 
done with replacement; this leads to the fact that about 
1/3 of sample data is sampled more than once. Boot- 
strapping has been applied to estimate cor relation func- 
tion e rrors in cosmology before, first by iBarrow et al.l 
(|1984f) , but it is intuitively clear that the bootstrap sam- 
ples represent a different world, where 1/3 of the galaxies 
of the original sample coincide (are close pairs) . So it has 
been avoided lately, and other methods (block jackknife, 
for example) are being used. 

The latter is a step in the right direction. The reason 
why the direct approach to bootstrap correlations fails is 
that bootstrap estimates can be app lied only to smooth 
funct i ons of samp le means (see, e.g. lEfron fc Tibshiranl 
[1991 ILahiril [20031) . Only in this case can the estimates 
be proved to be consistent (approaching the population 
statistics in the large sample limit). Correlation func- 
tions do not belong to this class. 

Ho wever, there is an elegant way to solve the prob- 
lem (jLahiril [2003D . Let us take a simple case of a ran- 



dom process Y(xi) defined on an onedimensional grid. 
Then, for studying its correlation function at the lag 
k, create another, two-dimensional process Z2{xi) = 
iY{xi),Y{xij^k))- Averaging over the product of the 
components of Z2 gives us the covariance function. 

Generalization of this approach to 3-D point processes 
involves introducing a fine 3-D grid, and generating a 
random process Y[xi) on the grid by assigning to its 
vertices values or 1 depending if the grid point (or cell) 
hosts a point of our point process. 

As above, let us create a new multidimensional process 
Zm{xi), a collection of m values oiY{-) at the grid points 
at a fixed distance interval (bin) from Xi . The number of 
points in the bin, if the original cell at Xi contains a point 
(the number of neighbours) can be defined as a simple 
function of this process, N{Zm{xi)). The sample mean 
of this function is the mean number of neighbours for a 
given distance bin, proportional to the function DD(s) in 
Eq. Ilj. This function can be bootstrapped. Notice that 
if wc consider all possible bins we get a full histogram of 
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neighbours of a sample point. Thus, the natural things 
to bootstrap are such pointwise histograms - histograms 
of all the neighbour distances from a given point, and, if 
needed, histograms of distances of random points from 
a given sample point. All points carry their histograms; 
generating a bootstrap realization of a correlation func- 
tion reduces to selecting a sample of histograms (with 
replacement), summing them to get the bootstrap real- 
ization, and applying, e.g., Eq. ([T]) to get the bootstrap 
realization of a sample correlation function. In order to 
be comparable to the original data, the number of his- 
tograms in the sum have to be the same as the number 
of points in the sample; as usual, about 1/3 of them oc- 
cur in the sum more than once. Random point pairs 
should not be bootstrapped, as these are merely a tool 
for performing volume integrals. 

This approach is logical and easy to apply. However, 
not even these histograms are independent, as required 
- the locations of sample points are correlated, and this 
has to be taken into account. The recipe for that is not 
to bo otstrap individual points, but data blocks (jLahiril 
|2003| ) . Such a procedure can be proved to give consistent 
estimates. The size of the block is determined by the 
correlation length. We shall show below that it is the 
case. 

An unusual (for cosmologists) feature of block boot- 
strap is that blocks can overlap; in fact, the estimates 
with overlapping blocks are more efficient than those 
with non-overlapping blocks. 

Numerica lly, this procedure is much cheaper to apply 
than that of lBarrow et al.l ()1984[ ) . If we have enough core 
memory, we calculate and store the pair distance proba- 
bility densities for every sample point, average them for 
block densities, and sum these for bootstrap realizations. 
As we do not need to recalculate distances, this boot- 
strap algorithm is fast. It is straightforward to modify 
this procedure for pairwise or pointwise weighting. 

4. APPLICATION TO POISSON-VORONOI 
PROCESSES 

We demonstrate and check the bootstrap machinery 
described above on a point field where the correlation 
function is known exactly - the Poisson-Voronoi process. 
As this process does not depend on the underlying ran- 
dom field, this is a clean and direct check of the recovery 
of the discre teness v ariance. 

Following Ivan de Weyga ert fc Ickd (|1989[ ). we start 
from a set of seed points distributed randomly in space 
(a Poisson point process). The set of all vertices of their 
Voronoi tessellation defines another (Poisson-Voronoi 
vertices) point process with a well known nume rically 
tract able two-point correlation function (Hcinric h et al.l 

[Hal). 

We used this process as a testbed for our correlation 
function estimates and the bootstrap recipes for its vari- 
ance and bias. For that, we generated 100 samples of 
Voronoi vertices within the DR7-LRG sample volume, 
with about the same mean density. We found correla- 
tion functions for all these realizations, estimated their 
bias and variance, and found the symmetric confidence 
regions. 

Fig. [5] shows the 98% and 80% confidence regions ob- 
tained from these realizations. Then we selected one of 
these realizations and looked how well the bootstrap pro- 
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Fig. 2. — Comparison of the sample variability and bootstrap pre- 
dictions for a Poisson-Voronoi process modelling the galaxy distri- 
bution in the DR7-LRG sample volume. The variability of a sample 
of 100 realizations of the process is shown by color, the 98% confi- 
dence region by red and the 80% confidence region by green. The 
block bootstrap estimates of variance are made on the basis of a 
single realization (its correlation function is shown by the yellow 
line), the predicted confidence regions are shown by blue lines (98% 
confidence) and by errorbars (80% confidence). The figure shows 
the overall fit of the observed and predicted variability, the inset 
shows that for a smaller region at a larger scale. 

cedure predicts the real bias and variance. The predicted 
variance was the closest to the true one for the block ra- 
dius R = 15 /i~^Mpc, which practically coincides with 
the correlation length (^(i?) « 1) for our process. The 
bootstrap confidence limits are shown in Fig. [5] by the 
blue lines (for 98%) and by errorbars (for 80%); the vari- 
ances are easier to compare in the inset. Biases were 
always considerably smaller than variances, both for the 
true samples and for the bootstrap samples. Based on 
these tests, we feel that the bootstrap procedure works 
well; we shall apply it to the galaxy correlation functions 
below. 

5. LARGE-SCALE GALAXY CORRELATION 
FUNCTIONS 

The redshift-space correlation functions for our sam- 
ples are summarized in Fig. |3l The correlation func- 
tions are found as described above, and their errors 
(biases, variances, confidence regions) are estimated by 
block bootstrap, with blocks as spheres with the ra- 
dius R = 12.0 Mpc for the DR7-LRG sample and 
R — 6.5 Mpc for the 2dFVL sample (the respective 
correlation lengths). We show in the bottom panel also 
the correlation function for a subsample of red galaxies 
from the 2dFVL (chosen according to the spectral type of 
the galaxy, see Madgwick ct al. 2002), with 17,252 galax- 
ies. The block size was chosen the same as for the main 
sample. 

The top panel shows the correlation function for the 
DR7-LRG sample. We see that the function ^(s) remains 
positive in the whole range shown(upto200ft."iMpc). It 
also shows confidently the 105 h~^Mpc maximum. The 
inset amplifies the large-distance behaviour of the LRG 
correlation function. We show the quantity s^^ (s) for the 
DR7-LRG and compar e it wi t h the correlation function 
estimates bv .Eisenstein et al.l (|2005l) . We show here only 
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Fig. 3. — The redshift space correlation function for the differ- 
ent samples. The top main panel shows the correlation function for 
the DR7-LRG sample (the 98% and 80% confidence regions in dark 
red and pale blue and the function as a green line). The acous- 
tic bump is clearly seen at about 105 /i~^Mpc. The inset shows 
the large-scale behaviour of the correlation function (multiplied by 
s^). The function is shown by the green line, its 80% confidence 
region in pale blue, and it is compa red with the original BAO 
discovery data llEisenstein et ani2005l ) with errorbars, labelled by 
Eis05. Note that the the acoustic peak is wider in the new sample. 
The bottom panel shows the redshift-space correlation function for 
the volume-limited samples: for the 2dFVL sample (the 80% confi- 
dence region in pink, the function in light green), for the red galaxy 
subsample of the 2dFVL (the red line with errorbars covering the 
80% confidence region), and for the DR7-LRG-VL (the 80% con- 
fidence region in pale blue, the function in dark green). The inset 
shows these correlation functions at short scales. Numerical values 
of the correlation functions plotted here can be downloaded from 
http : //www.uv . es/martinez. 

the 80% symmetric confidence region that is comparable 
with the ±cr limits of the previous estimate. We see that 
the new data shows the BAO peak with a much better 
confidence and also reveals that the peak is much wider 
than the peak found in D R3-LRG. A similar tr end was 
also seen in the analysis bv lCabre fc Gaztanagal (j2009). 

The bottom panel shows the 80% confidence level re- 
gion for the 2dFVL (by pink), and shows the 80% con- 
fidence limits for the red galaxies from the 2dFVL by 
errorbars. In contrast to the DR7-LRG, the correlation 
function for the 2dFVL sample crosses zero at about 
55 /i~^Mpc reaching a local minimum with ^(s) < 
at about 65 h~^Mpc. At larger scales, this correlation 



function tends also towards the 105 h~^yipc maximum. 
This maximum is present, even more prominently, in the 
correlation function of the red galaxy subsample. It is 
precisely at this scale where the LRG s ample shows the 
acoustic baryonic peak first noted by lEisenstein et al.l 
(120051 . 

The inset shows the behaviour of the correlation 
function of the volume-limited samples, DR7-LRG VL, 
2dFVL, and 2dFVL red, at short scales, characterised 
by a power-law regime with a downturn at the smallest 
scales (see, e.g., Ross et al. 2007). The difference in the 
correlation function amplitudes is a clear fingerprint of 
the luminosity-color segregation. The LRG galaxies are 
very luminous red objects lying in more dense environ- 
ments and displaying enhanced clustering. 

6. CONCLUSIONS 

The results show, first, that the baryon acoustic peak 
is a stable feature in the large-distance galaxy correla- 
tion function. We have presented the analysis of the 
LRG sample drawn from the last data release (DR7) of 
the SDSS. This survey is much larger and more compact 
than the data used before, so edge-correction effects are 
smaller now and results are therefore more robust. While 
the peak is most prominent in the distribution of lumi- 
nous red galaxies (because of the large spatial extent of 
the sample), it can be seen, if we know what to look 
for, in smaller galaxy samples. We have reported the 
first detection of the baryon acoustic peak in a volume- 
limited sample of very luminous galaxies drawn from the 
2dFGRS. 

Second, the baryon peak is much wider than found 
before and than expected; the distortions responsible for 
that demand careful analysis. 

Third, the minimum in the large-distance correlation 
functions of some samples demands explanation - is it 
really the signature of voids? 

And, fourth, a new internal method to estimate the er- 
rors of the correlation function has been introduced. The 
method is based on a generalization of the application of 
the bootstrap resampling techniques to smooth functions 
of sample means. 

Acknowledgements 

We thank Darren Croton for the 2dfGRS samples and 
the mask data, Rien van de Weygaert and Dietrich 
Stoyan for their advise regarding the Voronoi-Poisson 
process, Robert Lupton for useful discussions and Jun 
Pan and Jan Hamann for insightful comments on a pre- 
liminary draft of this letter. 

This work has been supported by the University of 
Valencia through a visiting professorship for Enn Saar 
and by the Spanish MEG projects AYA2006-14056 and 
GSD2007-00060, including FEDER contributions. PAM 
acknowledges support from the Spanish MEG through 
a FPU grant. ES and ET acknowledge support by the 
Estonian Science Foundation, grants No. 6104, 6106, 
7146 and by the Estonian Ministry for Education and 
Science, grant SF0060067s08. 

Funding for the SDSS and SDSS-II has been pro- 
vided by the Alfred P. Sloan Foundation, the Partic- 
ipating Institutions, the National Science Foundation, 
the U.S. Department of Energy, the National Aeronau- 
tics and Space Administration, the Japanese Monbuka- 



6 



Martinez et al. 



gakusho, the Max Planck Society, and the Higher Educa- tion Funding Council for England. The SDSS Web Site 

is http://www.sdss.org/. 

REFERENCES 



Abazajian, K. et al. 2008, arXiv:astro-ph/0812.0649 
Barrow, J. D., Bhavsar, S. P., & Sonoda, D. H. 1984, MNRAS, 
210, 19P 

Cabrc, A., & Gaztanaga, E. 2009, MNRAS, 393, 1183 

Cohn, J. D. 2006, New Astronomy, 11, 226 

Cole, S. et al. 2005, MNRAS, 362, 505 

CoUess, M. et al. 2001, MNRAS, 328, 1039 

Croton, D. J. et al. 2004, MNRAS, 352, 828 

Efron, B., & Tibshirani, R. J. 1993, An introduction to the 

bootstrap (London: Chapman and Hall) 
Eiscnstcin, D. J. ct al. 2001, AJ, 122, 2267 
Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 
Eisenstein, D. J. et al. 2005, ApJ, 633, 560 



Heinrich, L., Korncr, R., Mehlhorn, N., & Muche, L. 1998, 

Statistics, 31, 235 
Hiitsi, G. 2006, A&A, 449, 891 

Lahiri, S. N. 2003, Resampling Methods for Dependent Data 

(New York; Springer) 
Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 
Madgwick, D. S. et al. 2002, MNRAS, 333, 133 
Norberg, P., Baugh, C. M., Gaztanaga. E., & Croton, D. J. 2008, 

arXiv:astro-ph/0810. 1885 
Percival, W. J. et al. 2007, ApJ, 657, 51 
Ross, N. P. et al. 2007, MNRAS, 381, 573 

Saar, E., Martinez, V. J., Starck, J.-L., & Donoho, D. L. 2007, 

MNRAS, 374, 1030 
van de Weygaert, R., & Icke, V. 1989, A&A, 213, 1 
Zehavi, I. et al. 2005, ApJ, 621, 22 



